Energy spectrum of a 2D Dirac oscillator in the presence of a constant magnetic field 
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("^ , In this paper we obtain exact solutions of a 2D relativistic Dirac oscillator in the presence of a 

CNJ ■ constant magnetic field. We compute the energy spectrum and discuss its dependence on the spin 

,—1 ' and magnetic field strength. 
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^; I. INTRODUCTION 

<~i \ During the last years the study of electron in two-dimensional systems has become a subject of active research. 

A. . The rapid advances in microfabrication technology has made possible to confine laterally two dimensional electron 
q^ ' systems. These quantum confined few-electron systems are often referred to as artificial atoms where the potential of 
U , the nucleus is replaced by an effective potential. A parabolic potential V — 2^r'^ is often used as a realistic and at 
• ' the same time computationally convenient approximation. Despite its simplicity, parabolic quantum dots appear to 
d', be a good approximation to complicated quantum dot structures. Q 

^ • The solution of the two-dimensional Schrodinger equation for a single electron in a homogeneous magnetic field B 

has been known after the publication of the papers by Fock B and Darwin in the 20's. The inclusion of a parabolic 
'^ ■ potential to this problem leads to equations that can be solved without significative complications. 
5 ' Despite the large body of papers discussing the confinement of electrons in low dimensional structures, almost all 

of them tackle the problem in the framework of the Schrodinger equation. Among the difficulties when we consider 
relativistic Dirac electrons with parabolic confinement we have that we do not have at hand analytic solutions in 
terms of special functions and we need to use semi-analytical and numerical methods in order to compute the energy 
spectrum. 

The relativistic extension of a parabolic confining potential can be done with the of the Dirac oscillator [Q-g|. 
^vj ■ Among the advantages of this approach over the inclusion of Lorentz scalar potential ^mw^r^, as confining potential, 
ly-s [ we have the first one leads to a non-relativistic quadratic Hamiltonian in both coordinates and momenta and therefore 
f — its solution can also be expressed in terms of Laguerre or Hermite polynomials. The stability of the Dirac sea for the 
^^ Dirac oscillator Hamiltonian Q] also shows that this model does not present difficulties related to the Klein's Paradox. 
^~~l , In the present article we solve the two-dimensional Dirac oscillator in the presence of a constant magnetic field B 

^^ • perpendicular to the plane where the electron is confined to move. The article is structured as follows. In Sec. II, we 
solve the Dirac oscillator in polar coordinates. In Sec. Ill, we analyze the energy spectrum of the oscillator. In Sec. 
IV, we discuss discuss our results and the behavior of the energy levels in the non relativistic limit. Along the paper 
we adopt the natural units fi = I, c = 1. 

a 

O , II- SOLUTION OF THE DIRAC OSCILLATOR 

o 

The Dirac oscillator proposed by Moshinsky and Szczepaniak y includes a type of interaction in the Dirac equation 
which, besides the momentum, is also linear in the coordinates. The Dirac oscillator reduces, in the non-relativistic 
limit, to a harmonic oscillator with a very strong spin-orbit coupling term. Namely, the correction to the free Dirac 
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has the form: 

p ^ p — imujPr. (2) 

After substituting (0) into (IT]) we get an Hermitian operator linear in both p and r. Recently, the Dirac oscillator 
has been studied in spherical coordinates and its energy spectrum and the corresponding eigenfunctions have been 
obtained [^ . Since we are interested in studying the Dirac oscillator in a two-dimensional space, a suitable system of 
coordinates for writing the harmonic interaction are the polar p and i9 coordinates. In polar coordinates (i, r, "(9) the 
metric tensor gap takes the form 

9ap^diag{-l,l,r'^), (3) 

and the corresponding unitary vectors ip and e,? are 

Cp = cos 6i + sinOj, eg — —smOi + cos Oj (4) 

In this case the radial component of the modified linear momentum takes the form: pp — imujPp. The vector potential 
A associated with a constant magnetic field interaction is 

A^^e,. (5) 

The corresponding B can be written in polar coordinates as follows: 

B^Be, (6) 

Since we are interested in computing the energy spectrum of a two dimensional confined electron in the presence of a 
constant magnetic field, in the present article, we analyze the solution of the 2+1 Dirac oscillator in the presence of 
a constant magnetic field associated with the vector potential (||). 

One begins by writing the Dirac equation (|l|) in a given representation of the gamma matrices. Since we are dealing 
with two component spinors it is convenient to introduce the following representation in terms of the Pauli matrices 

/37i = CTi, /372 = 0-2, /3 = 0-3 (7) 

It is worth mentioning that Dirac matrices expressed in the Cartesian tetrad gauge take the form 

jP — ■J ■ Cp — j^ cos ^ + 7^ sin -d (8) 

7* = 7 • e^ = — 7I sin ?9 + 7^ cos ?9 (9) 

the similarity transformation S{{)) which reduces the Dirac matrices 7'' and 7 to 7^and 7^ is 

S{i})^cos--j^j^sin- (10) 

with 

S-\^)YS{^) = 7\ S~\^)j^S{^) = 7'- 

Using the matrix representation (|7|) and the similarity transformation ( |lO|) , we obtain that the Dirac oscillator in the 
presence of a constant magnetic field (0) has the form 
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where the spinor ^ is expressed in the (rotating) diagonal gauge. It is related to the Cartesian (fixed) spinor by means 
of the transformation S{p, 'd) |§] 

* = Vp^(i?)"1*c (14) 

with 

S{i))-^d = *c (15) 

The matrix transformation S^t}) can be written as follows 

S{^) = eM-^lW) = exp(-z^a3) (16) 

The Y^ factor has been introduced in ( |l4| ) in order to eliminate the term '^ > ^ '- > — . Noticing that ^(i^) satisfies 
the relation 

S{d + 2TT) = -S{^) (17) 

we obtain, 

'J'(?9 + 27r) = -*(t9) (18) 

so we have k^ — N + 1/2, where N is an integer number. The rotating Dirac spinor "^^ can be written in terms of 
^c as: 

We can label the quantum states of the Dirac oscillator in terms of eigenstates of the parity operator. In fact, analo- 
gously to the three dimensional Dirac oscillator S we have that the parity operator commutes with the Hamiltonian 

(0) 

[P,H]^^0 (20) 

where the parity operator acts on the Dirac spinor in the Cartesian gauge as follows: 

P*e(r) - /3*c(-r) = fT3*(-r) (21) 

In cylindrical coordinates, the reflection respecting to the origin r —>■ r is obtained via the rotation 

P^p d^^ + TT (22) 

Taking into account the relation (nj) between 'J and $c as well as ( p^ we readily obtain 

"^ ' ^ \ e^(fe^+i/2)'?^2(p) J ^ ' 

and consequently 

P*c(r) = (-l)'=^-i/2^c(r) (24) 

therefore, the parity of the energy eigenfunctions is given by (— l)'='S~i/2. Reminding that the total angular momentum 
J expressed in the Cartesian gauge takes the form 

^ = -(^ + ^), (25) 

and taking into account the representation (M) for the Dirac matrices we have that 



J*c(r) - {k^ - ^a=^)*c(r) (26) 

then, the eigenvalues of the total angular momentum are 

J = fctf T 2 (27) 

Using the representation (0) , the spinor equation (|ll|) can be written as system of two first order coupled differential 
equations, 

i{E - m)*i(p) = (A + ^ _ pmCj)^2{p), (28) 

dp p 

i{E + to)*2(/o) = (^ - ^ + pmcD)*i(p), (29) 

dp p 

where the frequency u) can be written in terms of the Larmor frequency as follows: 

eB 

u) = uj + sujl ~ ijJ ^ (30) 

2m 



substituting (p9|) into (ESI) and vice-versa we arrive at 



r ^2 (fc^)(fc^ T 1) , -.„, , ,^ 2-2 2 , .p2 2n' 

--^ 5 \- mojiZk^ ±1) ~ m uj p + (L — m) 



J; ) = (31) 



It is not difficult to see that the solution of the second order equation ( pl\ ) for '^i can be expressed in terms of 
associated Laguerre polynomials i^(a;) [||,|0[ as follows 

*i = ci cxp(-x/2)a;5(i/2+^)Lj;(a;) (32) 

where we have made the change of variables 

X = mujp^ (33) 

the parameter p is 

M = ±{kii - 1/2) (34) 

and the natural number n satisfies the relation 

^'~_"'' +(lTl)(2fc^-l) = 4n (35) 



Since the function ^i must be regular at the origin, we obtain that the sign of p in ( p2D is determined by the sign of 
k^. In fact, for fc^ > we have that ^i reads 

*i = ci exp(-a;/2)a;'=''/2L^''~^/2(a;) (36) 

substituting (|3g) into (^ we arrive at 

*2 = 2^c^^^^^eM-x|2)x^'-+'y^Ll^_\"\^) (37) 

where ci is an arbitrary constant. 

Analogously, we obtain that the regular solutions for fc^ < are, 

*2 = C2 cxp(-a;/2)a;-'=''/22,-fc-»-i/2(2;) (38) 

^^ ^ 2zc2^^^ exp(-x/2)a;(i-'^--)/2^y2-fc, (^) ^gg^ 

il, -|- m 

where Ci is a normalization constant The expression ( pq ) can be rewritten as follows 

£;2 - to2 ^ 4 [n - e(-fc^)(fc^ - 1/2)] (tow + — ) (40) 

where <d{x) is the Heaviside step function. 



III. STUDY OF THE ENERGY SPECTRUM 

From the relation ( [40|) it is clear that the energy spectrum of the 2+1 Dirac oscillator depends on the value of k^. 
Notice that for positive values of k^ the energy of the system has the form 

eB 

E"^ ^ m^ + An{muj + —) (41) 

For fc^ < we observe that, the states with (n ± /, fc^ — 1/2 ± I), where I is an integer number, have the same energy. 
In this direction there are some differences with the spherical Dirac oscillator H . Despite in both cases bound states 
are obtained, for the 2+1 Dirac oscillator the energy spectrum presents extra degeneracies only for negative values 
of k^. In order to get a deeper understanding of the dependence of the energy spectrum on the spin we can take 
the nonrelativistic limit of the Dirac equation (|l^). In order to do that, it is advisable to work with Eq. (|3l|). The 
Galilean limit is obtained by setting E = m -\- e, and considering e << m. Taking into account that the first two 
terms in Eq.(|3l|) are associated with the operator P^, we obtain in the nonrelativistic limit 

2m 2m 2 2m 2 

notice that Eq. (H3) corresponds to the Schrodinger Hamiltonian of a modified harmonic oscillator with an additional 
spin dependent term given by —ui{k^ ^ 2)' This contribution is proportional to the frequency of the oscillator plus 
the Larmor frequency. 

The energy spectrum of the 2+1 Dirac oscillator can be computed using the standard derivation followed by 
Moshinsky and Szczepaniak y,|ll|. In fact, using the Bjorken and Drell ||l^ representation for the Dirac matrices, 
the Dirac oscillator coupled to an electromagnetic field reads 

(S- m)*i = fT • [(p - eA) +imcjr]*2 (43) 

(£; + m)*2 = cr • [(p - eA) - imtjr]*! (44) 

taking into account that the vector potential for a constant magnetic field B is given by (0) and we are working in 
2+1 dimensions, we obtain 

(^2 - m^)* ^ {{p^ + rn^Lo'^p^) - 2mu) - 4mujL;,S:,) * (45) 

where 

L, = (rxp), 5, = ^ (46) 

finally, recalling the energy spectrum of a non relativistic harmonic oscillator we readily obtain that expression (|4 
takes the form 

(S^ - 771^) * = ( 2miD{2n + (fc^ ip -) + 1) - 2mu! - 2mu!{k^ T -) ) * (47) 

E^ -m^ = [4n + 2{-2ki, + l)e(-fci,)] mu! (48) 

result that coincides with expression ([40|). 

We can label the eigenstates of the parity operator in terms of the total angular momentum j. If we define 

e = -1 if parity is (-1)^ (49) 

e = 1 if parity is (— l)-*^^ 



in both cases we have k^ = J + f , therefore we have that the energy spectrum (MH) can be written in terms of the 
parity e as follows 

E'^ -m^ = [4n + 2(2fc^e + l)Q{eki))] mu) (50) 



From (48) we can see that the energy spectrum of the Dirac oscillator in the presence of a constant magnetic field 



presents a dependence on B that differs from that obtained with the help of the Pauli equation with a harmonic 



potential, [l^jljl ■ In this case the spin couples to the magnetic field via the term — |^ , and there is a contribution of 
B on the energy for positive as well as for negative spin orientations. 

From (Bq) it is straightforward to obtain the energy spectrum S states. For fc^ > we have 

eB 

E —m = Anrau) — Animio A — — ) (51) 

on the other hand, for k§ < 0, the energy spectrum satisfies the relation 

E^ -m^ = A{n+]--ki)){muj+^—) (52) 

from ( pl| ) and (m) we observe that the energy spectrum of the Dirac oscillator in the presence of a constant B depends 
on the spin orientation. For fc^ > the energy of S ground state satisfies E ~ m. Conversely, for fcij < the energy 
E satisfies the relation E^ = m^ + A{mu} + ^) and consequently depends on the magnetic field strength B. 

IV. DISCUSSION OF THE RESULTS 

The energy spectrum of a non relativistic two-dimensional electron confined by a parabolic potential in the presence 
a of static magnetic field B is |g,|| 

En' = (2n'+ I I I +\)hVt - -hujj (53) 

where n' = 0, 1, 2.., and / = 0, =Fl, =f2, .... The renormalized oscillator frequency O reads 



"=V^' + f (54) 

The expression ( |53| ) does not contain the spin field coupling term appearing in the Pauli equation: 

g*tiBB-S^±^B (55) 

4m 

that, for the value g* = 2, splits the energy levels by a factor ±Wc 

In order to better understand the energy spectrum (pl|), (p2), we proceed to analyze the asymptotic limits of the 
non minimal coupling introduced in Eq.(^). 

The weak oscillator limit can be obtained after taking the limit lu ^ 0. The energy spectrum (plh, ( p2[ ) reduces to 

E^ = m^ + 2[n~Q{-ki,){k^-l/2)]eB (56) 

which are the relativistic Landau levels obtained when one solves the Dirac equation. The non relativistic limit can 
also be obtained in a straightforward way from ( pq ) giving as result the expression ([5^). The energy spectrum in weak 
magnetic field limit presents a behavior that is not observed in the non relativistic oscillator case. When we eliminate 
the magnetic field strength B the energy spectrum remains degenenerate but depends on the value of k^ therefore, in 
this limit the Dirac oscillator does not seem to be a suitable approximation to the problem of an electron confined in 
a quantum dot [|l],^ . 

V. CONCLUDING REMARKS 

The presence of a constant magnetic field B (j^) does not break up the symmetry of the vacuum and therefore 
we do not have mixing of positive and negative energy solutions. The confining Dirac oscillator potential gives as a 
result a spin orientation dependent energy spectrum for the electron in a the presence of a constant magnetic field. 
There are extra degeneracies of the quantum states with the same energy when fc^ < O.The non relativistic limit of 
the model presented in this paper shows that B modifies the oscillator frequency uj via the Larmor term ujl = ■^■■ 
The nature of of the coupling associated with the Dirac oscillator does not not permit one to use this model in the 
study of an electron confined by a quantum dot. 
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